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Abstract We study 7530 sunspot umbrae and pores measured by the Hinode 
Spectropolarimeter (SP) between November 2006 and November 2012. We pri¬ 
marily seek confirmation of the long term secular decrease in the mean magnetic 
field strength of sunspot umbrae found by Penn and Livingston (2011, lAU 
Symp. 273,126) between 1998 and 2011. The excellent SP photometric properties 
and full vector magnetic field determinations from full-Stokes Milne-Eddington 
inversions are used to address the interrelated properties of the magnetic field 
strength and brightness temperature for all umbral cores. We find non-linear 
relationships between magnetic field strength and umbral temperature (and con¬ 
tinuum contrast), as well as between umbral radius and magnetic field strength. 
Using disambiguated vector data, we find that the azimuths measured in the 
umbral cores reflect an organization weakly influenced by Joy’s law. The large 
selection of umbrae displays a log-normal size spectrum similar to earlier solar cy¬ 
cles. Influenced by the amplitude of the solar cycle and the nonlinear relationship 
between umbral size and core magnetic field strength, the distribution of core 
magnetic field strengths, fit most effectively with a skew-normal distribution, 
shows a weak solar cycle dependence. Yet, the mean magnetic field strength 
does not show a significant long term trend. 
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1. Introduction 

Despite their often complex configurations and dynamic evolution, simplistic 
properties of sunspots when viewed in collection set important constraints on 
the activity cycles of our Sun, as well as other stars. The record of the waxing 
and waning number of individual sunspots and sunspot groups forms one of the 
longest astrophysical studies to date. Sunspots, of course, are not individual. 
Subject to Gauss’s law for magnetism, sunspots form within bipolar active 
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regions that emerge from below the photosphere to establish paired groups 
of opposite polarity spots and/or plage. This process is organized according 
to Hale’s polarity law, Joy’s law for the tilt of bipolar active regions, and an 
equator-ward drift of a northern and southern band of active latitudes (see. 


e.g. Stix (2004)). Each of these solar cycle characteristics inform us about the 


mechanisms generating stellar activity; yet, much is unknown, including what 
leads to long-term declines in activity such as the Maunder Minimum between 


about 1645 to 1715 (Eddy, 1976). 


Recently, there has been renewed interest in the average magnetic and thermal 
properties of sunspots over the solar cycle. Measurements from the McMath- 
Pierce Telescope of the magnetic field strength and continuum brightness within 
the darkest portion of sunspot umbrae {i.e. umbral cores) exhibit on average a 
long-term decline in field strength and a long-term increase in brightness between 


1998 and 2011 ( 

Penn and Livingston, 2006 

Penn and Livingston, 2011 

Liv- 

ingston, Penn, and Svalgaard, 2012 

). These data examine the Zeeman-splitting 


of the Stokes-I a components of the Fe I normal Zeeman triplet at 1564.8 nm 
(Lande factor, g = 3). In contrast, a study of archived magnetic field measure¬ 
ments gleaned from polarized states of visible spectral lines from the former 
USSR indicates that the maximum field strength in sunspots rise and fall during 
the solar cycle and do not show a long-term decline between 1957 and 2011 


(Pevtsov et ai, 2011). This cyclic behavior is supported by the line-of-sight 


(LOS) umbral magnetic flux measurements from Watson, Fletcher, and Marshall 


(2011). The primary difference between these studies and that of Livingston, 


Penn, and Svalgaard (2012) is that only the largest spots on any given observing 


day are included in the temporal averages of the umbral maximum magnetic field 


strength by [Pevtsov et al. 

(2011 

) and 

Watson, Fletcher, and Marshall 

(2011) 

while 

Livingston, Penn, and Svalgaard 

(2012 

) include all observed spots. 


Determining whether the average thermal and magnetic properties of sunspots 
vary in time intricately depends upon the accurate characterization of the dis¬ 
tributed magnetic field, thermal, and size parameters of sunspots. The maximum 
umbral field strength and minimum intensity are closely related with the size 


of the sunspot umbra (Kopp and Rabin, 1992 Mathew et ai, 2007 Schad 


and Penn, 2010), while the maximum magnetic field strength and minimum 


intensity of sunspot umbrae are well correlated with each other (Norton and 


Gilman, 2004 Schad and Penn, 2010). Nagovitsyn, Pevtsov, and Livingston 


(2012) suggest that, both the solar cycle variation of the maximum magnetic 


field of the largest sunspots during a solar cycle inferred by [Pevtsov et al. (2011) 
and the long-term decline in the average sunspot magnetic field reported by 


Penn and Livingston (2006), can be explained by temporal variations in the 


relative distribution of small and large sunspots, as seen in full disk synoptic 
images recorded at the Kislovodsk Mountain Astronomical Station of Pulkovo 
Observatory. This is in disagreement with studies of the relative size distribution 
of umbrae measured between 1917 and 1982 using Mount Wilson whitelight 


images (Bogdan et ai, 1988), which exhibits a log-normal distribution that 


varies only in amplitude during the solar cycle. The invariant log-normal size 


distribution (Bogdan et ai, 1988) is further supported by measurements from the 
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Greenwich Observatory (Baumann and Solanki, 20051 and the Kitt Peak Vacuum 
Telescope’s Spectromagnetogragh (KPVT/SPM) (Schad and Penn, 2010). 

Schad and Penn (2010) and Mathew et al. (2007) discuss the role of the 


size dependence of the umbral intensity measurements to argue against a solar 
cycle dependence of the average umbral core intensity. [Schad and Penn] also 
show, using LOS magnetic flux measurements, that in an automated selection 
of sunspot umbrae between 1993 and 2004, the average magnetic held strength 
of sunspot umbrae does not exhibit any signihcant temporal changes, nor do 
the relationships between umbral core intensity, size, and maximum magnetic 
held strength. The role of the size dependence is not addressed by [Penn and 


Livingston (2006), who used the McMath-Pierce Telescope data set, as umbral 


sizes are not measured. Although the authors assume that the measurements are 
not biased by selection ehects, this is not explicitly demonstrated. 


A great disadvantage of the studies by Schad and Penn (2010) and Watson, 


Fletcher, and Marshall (2011) is, in contrast, the primary advantage of the Penn 


and Livingston (2006) data set. The Fe I 1564.8 nm intensity spectra is sensitive 


to the total magnetic held strength within the sunspot umbrae, since its Zeeman 
(T components are well resolved in most umbrae. Umbral magnetic held measure¬ 
ments from KPVT/SPM and the Michelson Doppler Imager on board the Solar 
and Heliospheric Observatory (SOHO/MDI) only measure the LOS magnetic 
hux, and are subject to other errors inherent to the intensity-difference and 
center-of-gravity magnetograph techniques. Furthermore, to estimate the total 
magnetic held in umbral cores, one must make the assumption of a unity hll 
fraction and that the umbral magnetic held is vertical {i.e. radial) within dark 
umbral cores, which until now has not been tested for a large collection of umbrae 
(see Section]^. 

Modern full-Stokes spectropolarimeters infer the true umbral magnetic held 
strength via full Stokes analysis of spectral lines. Although no synoptic full- 
disk vector measurements are available for the full time period of Penn and 
Livingston’s ? measurements, the six years (November 2006 - present) of high 


resolution full-Stokes measurements from the Spectropolarimeter (SP) (Tsuneta 


et al, 20081 

Lites et al., 2013 

Lites and Ichimoto, 2013) onboard the Hinode 

spacecraft ( 

Kosugi et al, 2007 

offers a unique chance to investigate the collec- 


tive thermal and magnetic held vector properties of sunspot umbrae. Below we 
discuss the full vector structure of dark umbral cores measured by Hinode/SP, 
presenting in detail the interrelated magnetic held, thermal, and size parameters 
of sunspot umbrae (see Sections and . We hnd an agreement between our 


measured distributions of umbral sizes and the study of Bogdan et al. (|1988), 


and, thus, we argue that our selection adequately samples the solar cycle (see 
Section 1^. The core magnetic held strengths are skew-normally distributed with 
weakly signihcant variations over the solar cycle; however, the average magnetic 
held strength does not display any long term decrease (see Section!^. 


2. Observations 

The SP onboard the Hinode spacecraft regularly achieves a polarimetric accu¬ 
racy of 10“^ of the incoming intensity when measuring the full Stokes polarized 
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spectra of the Fe I 630.15 and 630.25 nm spectral lines with a spectral sampling 
of 2.15 pm (Kosugi et al, 2007 Tsuneta et ai, 2008 Ichimoto et al.^ 2008). The 
SP spectrograph slit has a projected length of 164" and spatially samples the 
solar surface every 0.16" both along the slit and in the stepping direction, which 
is restricted to 328” (Tsuneta et al., 2008). As a proposal driven instrument, 
the mode of operation and target selection for Hinode varies. All of the spatial 
SP maps analyzed here have at minimum an effective resolution of 0.64", cor¬ 
responding to a spatial sampling of 0.32". The SP is a stable platform. Seeing 
fluctuations are non-existent and platform vibrations do not degrade the quality 
of the data. 


2.1. Data Selection and Azimuth Disambiguation 


We analyze 628 Hinode/SV area scans, manually selected from the Level 2 data 
archive between November 2006 and November 2012. The SP often scans the 
same region many times during its progression across the solar disk. Further¬ 
more, larger active-region targets are typically preferred over smaller sunspots 
within the Hinode target selection process, which might introduce a selection 
bias if all scans are included. To avoid this, we do not examine every area scan. 
Instead, on each UT calendar day, we select an individual scan for each unique 
targeted region, and attempt to maximize the time interval between scans taken 
on adjacent UT days. Scans that have a large field-of-view are preferred in our 
selection over scans of a reduced field-of-view, which typically were acquired 
in a higher resolution mode. We also ensure that the scans have full coverage 
over each individual sunspot umbra. Of the selected scans, 72 had a high East- 
West (x) (North-South (y)) spatial sampling of 0.148"(0.160"), 553 had an x(y) 
sampling of 0.297"(0.320"), and 3 had a mixed x(y) sampling of 0.148"(0.320"). 

The Level 2 data processing includes inversions of the measured Fe I polarized 
spectra at each spatial pixel using the MERLIN Milne-Eddington inversion code 
(Lites, 2009). MERLIN fits the intensity spectrum with a weighted combination 


of a single magnetic atmospheric component and a stray-light component using 
the normal definition for the fill fraction. The magnetic component is charac¬ 
terized by a single magnetic field strength, and can be considered an average 
value over the range of spectral-line formation in the atmosphere. The stray-light 
profile is determined from the average intensity spectrum from regions within 
the scan with low polarization. For sunspot umbrae, the magnetic fill fraction 
is typically high (> 0.8). To study the inclination of each umbral core magnetic 
field in the local reference frame of the sunspot, we resolve the 180° azimuth 
ambiguity of the transverse Zeeman effect within each sunspot umbra for the 
159 target scans located near the central meridian (±10°) using the AZAM 
disambiguation tool ( Lites et al, 1995] ) prior to transforming the LOS vector into 
the local reference frame. Due to the divergence (convergence) of the sunspot 
magnetic fields for positive (negative) polarity umbrae, the disambiguation of 
the umbral fields is a simple, straightforward process. 


2.2. Automated Threshold-based Sunspot Detection 


A total of 7530 dark features are found in our Hinode data set using an au¬ 
tomated method based on intensity thresholding of the continuum intensity 
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Figure 1. {top left) An example of our automated selection of sunspots and pores for one 
scan from the Hinode Spectropolarimeter Level 2 data archive. All structures with a normalized 
intensity below 0.575 are identified as a sunspot feature. The total magnetic field {bottom left) 
is recorded in these Level 2 data from a Milne-Eddington inversion, which includes weighting 
by a non-magnetic filling factor determined directly from the polarized spectral fits. The yearly 
distribution of all 7530 selected spots adequately samples the latitudinal dependence of sunspot 
emergence through the solar cycle as shown in the butterfly diagram of the core magnetic field 
strength {bottom right), while the total number distribution tracks the solar cycle variation in 
the SIDC international sunspot number {top right). 


maps for each scan, as in Schad and Penn (2010). These dark features consist of 
sunspot umbrae as well as pores without penumbra, which we collectively refer 
to as simply umbrae. The continuum intensity of each scan is first corrected 
for the center-to-limb dependence of the continuum intensity at 630 nm using 
the wavelength-dependent values given by Pierce (2000). The intensity is then 
normalized over the entire scan by the 75th percentile value of the intensity; 
that is, we divide by the value of the center-to-limb corrected intensity that 
separates the upper 25% of the intensity values from the lower 75% values. 
We find this a better automated method for normalization due to the varying 
amount of plage. Regions identified as umbrae are required to have a normalized 
intensity value below 0.575. Dilation and erosion image processing techniques 
are employed to connect neighboring regions separated by a small gap. Every 
selected feature must also contain at least 4 observed spatial pixels. An example 
of the automated selection is presented in Figure 

The slit-scanning mirror responsible for scanning the solar image across the 
Hinode/SP spectrograph slit does not translate linearly (Centeno et ai, 2009). 
Instead, it periodically deviates from a linear translation with an amplitude 
of about one arcsecond. While this small error is not expected to introduce 
large errors in the location of the darkest portion of the umbra from which 
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we glean parameters, it does influence the determined umbral sizes. We correct 
the non-linear scanning using the calibrated map labeled by the ‘X_coordinate’ 
keyword within the Level 2 data. Spatial masks recording the identified umbrae 
are destretched to the median scanning step prior to measuring umbral sizes. 

In total, 1117 umbrae are found in the high resolution maps, 6408 in the lower 
resolution maps, and 5 in the mixed resolution maps. We record the location, 
normalized intensity (or contrast), magnetic field strength, LOS inclination, and 
LOS azimuth for both the darkest pixel in each umbra as well as the pixel 
with the greatest magnetic field strength. We find a good correlation between 
the location of the largest field strength and darkest intensity. Thus, we limit 
our study only to the pixels with the darkest continuum intensity within each 


umbra in accordance with previous studies (. 

^enn and Livingston, 2006 

Penn 

and MacDonald, 2007 

Schad and Penn, 201( 

)). The total area of each umbra 


is determined by the number of spatial pixels it covers and the foreshortened- 
corrected effective area of each pixel in physical units. The radius of each umbra 
is then calculated assuming the umbra is circular. 

As shown in Figure [^, the total number of spots found in our data selection 
varies with the solar cycle in accordance with the Solar Influences Data Analysis 


Center (SIDC) sunspot number (SIDC-team, 2012). The magnetic butterfly dia¬ 


gram in Figure gives the core magnetic field strength and latitudinal location 
of each umbra as a function of the solar cycle. Low latitude umbrae from Solar 
Cycle 23 are measured between November 2006 and May 2008. Umbrae from 
Solar Cycle 24, starting in late 2008, exhibit the familiar equator-ward drift in 
latitudinal locations as a function of solar cycle. We conclude from this figure 
that our selection of spots adequately samples umbrae from the full latitudinal 
extent of sunspot formation during the solar cycle. 


2.3. Stray-Light Considerations 


We correct the determined umbral core intensities for instrumental stray light 
by accounting for the fill fraction values reported directly by the MERLIN 


inversions. According to Lites (2009), during the inversions the Stokes spectra 


are fitted with a single magnetic component and a field-free stray-light profile 
according to: 


/sy„(A,/,x) = //mag(A,x) -k (1 - f)Istrav{A) (1) 

where Isyn denotes the synthetic Stokes profiles fitted to the observed profiles. 
These calculated profiles result from the weighted contribution of a magne¬ 
tized atmospheric component, Imag, and a scattered light profile, IstrayW The 
properties of the magnetized atmosphere are described by the standard set of 
parameters, x, including the Milne-Eddington factors, thermal line width and 
damping, and a single magnetic field vector. The weighting factor, f, is known as 
the magnetic fill fraction or, alternatively, (1 — /) denotes the stray-light fraction. 
I stray {^) Is determined directly from the data by averaging the profiles within 
spatial points of low total polarization. Then, the returned magnetic field already 
has taken into account the role of stray light, where as our determination of the 
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intensity has not. We correct the umbral contrast measurements, as in Mathew 
et al. (20041, by dehning, in accordance with Equation Q , the corrected umbral 
contrast as: 


DC _ 

^spot 


^Ibs + (/ ~ ^)R‘stray 

f 

Kbs + (/-!) 

/ 


( 2 ) 


where is the umbral intensity contrast measured in the continuum, deter¬ 
mined directly from the HINODE/SP intensity maps that have been flattened 
and normalized. We make the approximation that the light scattered into the 
umbra, originates from the quiet Sun, such that Rgtray = 1- 


3. A Continuum of Thermal and Magnetic Field Properties 
3.1. Relationship of Umbral Contrast to Magnetic Field Strength 


King (1934) first discovered the strong correlation between the peak magnetic 


field strength and the photometric intensity of a collection of sunspots; further¬ 
more, he discussed the wavelength dependence of this correlation. Since then, 
the empirical relationship between the emergent intensity and the magnetic field 
strength has since been studied by a number of authors both within individual 


sunspots (Martinez Pillet and Vazquez, 1993 Solanki, Walther, and Livingston, 


1993 Leonard and Choudhary, 2008 

Jaeggli, Lin, and Uitenbroek, 2012 

) and 

amongst a collection of sunspots ( 

Kopp and Rabin, 1992 

Penn et al, 2003 

Nor- 

ton and Gilman, 2004[ |Schad and Penn, 2010 

Rezaei, Beck, and Schmidt, 


2012). Collections of sunspot umbrae observed by the SOHO/MDI (Norton 


and Gilman, 2004) and the KPVT/SPM (Schad and Penn, 2010) exhibit power 
law relationships between the minimum umbral intensity and the corresponding 
magnetic field strength. Yet, these studies are hindered by low number statistics 


(esp. Norton and Gilman (2004)), and scatter-inducing instrumental effects {e.g. 


stray-light and magnetograph “saturation” effects). 

Taking advantage of the excellent photometric and spectropolarimetric prop¬ 
erties of the SP, we re-investigate the relationship between intensity and mag¬ 
netic field strength for our collection of 7530 identihed umbral cores (see Fig¬ 
ure]^. The relationship between the magnetic field strength and the continuum 
contrast is shown both for all of the umbrae (top panels in Figure /r > 0.45, 
where p. is the cosine of the umbra’s heliocentric angle) and for the spots close 
to disk center (bottom panels in Figure]^ p. > 0.95). Evidenced by the lack of 
linearity in the right log-log scaled plots, a single power law relationship fails to 
accurately characterize these data. Instead, we empirically describe the trends 
via three linear fits for the ranges of R indicated in the figure. Furthermore, 
the scatter in these data is greatly reduced in comparison with previous studies. 
Since the SP resolves individual umbral dots, we expect our ability to locate the 
darkest portion of the umbra to be improved, which, in addition to a reduced 
stray light, may explain a lower scatter. 


SOLA: schad_revised2.tex; 22 May 2015; 0:16; p. 7 



















































Schad 


Temp (t63o = 1 )[K] 

3500 4500 5000 5500 


4000 


3000 

<0 

3 

ca 

S 2000 

m 


1000 


0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Umbral Core Intensity, Ic/Ic,qs 



Temp (t83o = 1 )[K] 

3500 4000 4500 5500 6182 



Umbral Core Intensity, Ic/lc.os 


Temp (t63o = 1 )[K] 

3500 4500 5000 5500 


4000 


3000 

¥ 

3 

to 

55 . 2000 

CQ 


1000 


0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Umbral Core Intensity, Ic/Ic.qs 



Temp {Te3o = 1 )[K] 

3500 4000 4500 5500 6182 



Umbral Core Intensity, Ic^Icqs 


Figure 2. Magnetic field strength versus normalized intensity for dark umbral cores. The 
top panels show the relationship for all the identified spots, while the bottom panels show 
the relationship for umbrae measured near disk center. Left and right plots correspond to 
linear-linear and log-log representations of the same data. Linear fits (blue continuous lines) 
of the observed data are over-plotted, and quantified within the figure for three ranges of the 
normalized intensity (denoted by R). Note that the magnetic field strengths for the darkest 
umbrae (R< 0.06) are anomalously high due to the effect of molecular blends on the Stokes 
inversion process. The temperature scale given at the top of each panel refers to the umbral 
core temperature at unit continuum optical depth inferred using Planck’s law. 


For high-contrast spots, a pronounced, nearly-isothermal strengthening of the 
magnetic field amongst the darkest umbral cores is evident. However, while 
such behavior has been observed within individual spots using the infrared Fe 
I lines at 1565 nm (Jaeggli, Lin, and Uitenbroek, 20121, the behavior seen here 
is probably a result of the influence of TiO and MgH molecular lines on the 
MERLIN inversions. We extract the Hinode Stokes spectra from the 29 scans 
with an umbra exhibiting a normalized intensity below 0.1 and a magnetic field 
strength above 3250 G. This subsample contains only 12 unique umbrae. We 
perform independent Milne-Eddington inversions on the spectra of the darkest 
umbra. The results return inconsistent measures of the magnetic field. Though, 
in general, the inferred magnetic field strengths for these umbrae are reduced 
in relation to the Level 2 Hinode values. It is clear that these spectra must 
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Figure 3. Squared magnetic field strength versus umbral core brightness temperature for 
all the measured umbrae (/i > 0.45, left panel), and for those umbrae close to disk center 
> 0.95, right panel). The linear fits found in Figurej^are over-plotted, using Equation 
to convert normalized intensities into brightness temperatures. 


be analyzed with more appropriate atmospheric models incorporating molecular 


formation and the molecular Zeeman effect (see, e.g., Wenzel et al. (2010) and 


Berdyugina (2011)). 


3.2. Temperature - Magnetic Field Relationship 


A theoretical understanding of the intensity-magnetic field strength relation¬ 
ship requires a radiative magneto-hydrodynamic (MHD) approach akin to the 
studies of Rempel et al. (2009) and Borrero, Rempel, and Solanki (2010). Yet, 


a simple picture can be built from the static force balance consisting of the 
external gas pressure from the photosphere/upper convective zone, the internal 
gas pressure of the sunspot, and the horizontal component of the magnetic 


Lorentz force. Following the assumptions of Maltby (1977) and the notation 


of Solanki, Walther, and Livingston (1993), the horizontal magnetohydrostatic 


(MHS) pressure equation can be written in cgs units as: 


Po{z) = Pg{r,z) + 


BUr,z) + Fc{r, z) 
Stt 


where 




( 3 ) 


( 4 ) 


Pq refers to the external gas pressure, while Pg is the internal gas pressure of 
the sunspot, r is a radial coordinate measured from the sunspot center in a 
cylindrical geometry with z being the height coordinate. The Lorentz force is 
decomposed into the two terms corresponding to a pressure term and the field 
curvature tension term, Fc. As in [Jaeggli, Lin, and Uitenbroek (2012), under the 
limited assumptions of a vertical magnetic held in a held-free, ideal background 
atmosphere of constant height and density. Equation ([^ reduces to a simple 
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proportionality between the squared magnetic field strength and the difference 
between the external and internal temperatures, 


i.e. 




OC Tqs - Ts 


(5) 


Accordingly, we display in Figure the relationship between the squared mag¬ 
netic field strength and the umbral core temperature determined from our col¬ 
lection of observed umbral cores. 

The umbral core temperature is determined from the normalized continuum 


intensity and the Planck function, as in Solanki, Walther, and Livingston (1993), 
using the equation 


DC _ 

■^spot 


^ spot 


I: 


QS 


^hc/XkgTo _ ^ 
ghc/XksT _ 2 ^ ’ 


( 6 ) 


where h, c, A, ks take their customary meanings. Igpot and Igg are the intensities 
determined in the spectral continuum within the sunspot and the local quiet 
Sun, respectively. T and Tg are the brightness temperatures of the spot and 
of the quiet photosphere. Tg is assigned for our 630 nm channel to 6182 K, in 
accordance with the measurements in Maltby et al. (19861. 


Just as in the studies of the temperature-magnetic field strength relationship 
within individual sunspots, dark umbral cores when viewed in collection do not 
comply with the simple picture of MHS balance given by Equation ([^. As shown 
in Figure while the temperature and squared magnetic field strength are 
clearly well correlated, the relation is not linear. Once again, the strengthening 
of the umbral field for umbrae with low temperatures produces a pronounced 
effect; although, this is probably an artifact induced by the molecular formation 
(see the discussion in the previous section). Furthermore, umbrae with higher 
temperatures exhibit a non-linear relation. As particularly apparent in the disk 
center sampling of umbrae, the slope of the vs. T relationship varies between 
intermediate and warmer temperatures, which is not found in other studies. It is 
unclear at this time what generates this difference; though, we note that we do 
not distinguish in our sample between umbrae with and without penumbra. The 
smaller, brighter features that we measure probably have no penumbra, while the 
larger, darker features do. The presence of a penumbra can significantly alter the 
force balance throughout a magnetic field concentration, which might influence 
the B^ vs. T relationship. Furthermore, the formation height of the continuum 
within sunspot umbrae is influenced via the well-known Wilson depression ef¬ 
fect by the magnetic field intensity. Consequently, the simple proportionality 
of Equation ^ no longer holds, leading to changes in the observed B'^ vs. 
T relationship. Measurements of the Wilson depression depth can be used to 
constrain the curvature force within an umbraa, but such measurements are 
difficult to achieve. 


4. Umbral Size Dependencies 

In Figure the magnetic field strength and normalized continuum intensity 
measured in the umbral cores are related to the overall umbral size, quantified 
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Figure 4. Dependence of the umbral core intensities (left panels) and the total magnetic 
field strengths (right panels) on the umbral radius (top panels) and umbral area in square 
megameters (bottom panels). Linear-log and log-log regression curves (blue continuous lines) 
provide empirical fits to the data points. The equations within the panels are in the same units 
as the left and bottom axes of each plot. As a reference, umbral areas are given in units of 
microhemispheres in the top axes of each panel. 


both as a radius (he. an “effective radius” assuming a circular umbra) and the 
total fore-shortened corrected area. Even though the Hinode/SP maps are of 
considerably higher resolution than the data used in previous determinations 
of the temperature-magnetic size dependence, considerable scatter still exists 
can be seen in the figure. The complex geometry of many umbral features is no 
doubt at the origin of this scatter. A regressive fit is shown for each relationship in 
Figure Q. We emphasize that umbral size (as well as intensity and magnetic field 
strength) depends on the definition of the umbra-penumbra boundary intensity 
and the wavelength of observation. Yet, it is clear from these observations that 
the relationship between the umbral size and the thermal-magnetic properties of 
the umbra is not linear, which has interesting consequences below (see Section 

§. 
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Figure 5. Histograms of the unsigned deviation of the umbral core magnetic field inclination 
from the local vertical axis. As expected, weaker umbral cores present a larger deviation range. 
TjN give the total number of umbrae included within each histogram. 


5. Deviation from the Vertical Field Assumption 
5.1. How vertical are magnetic fields in umbral cores? 


The angle between the magnetic field vector and the LOS direction is quantified 
as the magnetic field inclination angle (in the LOS reference frame). Zeeman 
diagnostics deterministically return the inclination angle via the analysis of the 
Stokes profiles. Unlike the azimuthal angle {i.e. the projection of the magnetic 
field vector on the plane perpendicular to the LOS direction), the inclination 
angle determination is unambiguous. Yet, this inclination angle is measured with 
respect to the LOS, not to the local solar reference frame. In order to determine 
the inclination in the solar reference frame, a transformation, which involves not 
only the inclination angle but also the azimuthal one, must be performed. Thus, 
the determination of the local inclination angle of the magnetic field is subject 
to our ability to resolve the azimuthal angle ambiguity in the LOS frame. 

For a quick estimation of the deviation of the umbral magnetic field vector 
orientation from the local solar vertical, one can calculate the deviation of the 
measured inclination angle in the LOS coordinate system from the expected 
inclination angle of the magnetic field in this same reference system, assuming 
a perfectly vertically directed field vector. This method evaluates the deviation 
of the core umbral magnetic fields from the vertical direction without the need 
for azimuth disambiguation. This quantity can be written as 


AO 


vertical 


arccos (fi) — arccos 


\B cos 0b\ 
B 


(7) 
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Figure 6. An illustration of the azimuthal angle reference frame used in this work, and Hale’s 
and Joy’s law for Solar Cycle 24 (December 2008 - present). The diagrams on the left show 
a diverging and converging transverse magnetic field extending from a positive and negative 
polarity region, respectively. In the local reference frame of each sunspot, an azimuthal angle 
of 0° is directed towards solar West. Hale’s Law, as depicted in the right diagram, predicts 
that negative polarity spots would lead in bipolar active regions in the northern hemisphere 
for Solar Cycle 24. 


where /r is the cosine of the angle between the LOS and the local solar vertical, 
B is the magnitude of the magnetic held, ds is the inclination of the magnetic 
held in the LOS reference frame that ranges from 0° to 180° in these Level 
2 Hinode inversions. The second term in Equation Q may be thought of as 
an “unsigned LOS magnetic held inclination angle”, whose range is between 0° 
and 90°. Histograms of this deviation angle are given in Figure]^ The one- 
sigma deviation from the vertical held approximation for strong umbral cores is 
^ 9°, conhrming earlier results of vertically directed helds in strong umbra. The 
vertical held approximation quickly vanishes for weaker umbrae. 

5.2. Joy’s Law and the Vector Magnetic Field of Umbral Cores 


Umbral core magnetic helds exhibit a range of deviation angles from the local 
vertical direction. We investigate if this deviation might be organized in any 
way. Karachik, Pevtsov, and Abramenko (20101 found that the peak magnitude 
of the average LOS magnetic hux for each polarity of a bipolar active region 
was measured on separate sides of the solar meridian, putting in evidence that 
opposite polarity magnetic helds within an active region may be tilted towards 
each other at the solar surface. We investigate this possibility using the infor¬ 
mation of the disambiguated magnetic held vectors of our identihed umbrae. 
The individual active regions are not investigated separately; rather, we test 
the key phenomenological laws describing the global organization of the solar 
magnetic held, i.e. Hale’s and Joy’s laws. Restricting the investigation only 
to spots observed after December 2008 within Solar Cycle 24, our sample is 
separated into four categories based on Hale’s polarity law for Solar Cycle 24, as 
shown in Figurej^ where spots are classihed by polarity and hemisphere. Umbral 
cores with local inclination angles less than (greater than) 90° are considered 
positive (negative) polarity spots. 
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Figure 7. Histograms of the local solar inclination angle of the umbral core magnetic field 
for umbrae classified according to Hale’s polarity law in Solar Cycle 24. The blue curve gives 
the raw histogram for 2° inclination bins, while the grey filled curves reports the probability 
density weighted by the relative area of each annulus of inclination. Gaussian functional fits 
to the raw histograms are performed, with the parameters given in each plot. 


Histograms of the local solar inclination angle of the umbral cores in the four 
categories are illustrated in Figure They are ordered according to the Cycle 
24 Hale’s polarity law with histograms of the northern hemisphere umbrae at 
the top. Umbrae of all sizes and strengths are included. Two distributions are 
included in the figure. The best representation of the number of umbrae with a 
given local inclination value is given by the regular histogram (blue curve) for 2° 
bins. The parameters describing the fit to these curves by Gaussian distribution 
functions are given in the figure. Note that these distributions have peaks away 
from 0° and 180°, which correspons to radially directed fields. This offset does 
not imply that umbral core magnetic helds are preferentially non-radial. Rather, 
to calculate a probability density function for the umbral core inclination angles, 
one must take into account the relative area of each inclination annulus in a 
spherical geometry. To do this, we calculate the density function: 


PDF{Ae) = 


1 ANe 
pAA/A' 


( 8 ) 
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Figure 8. Histograms of the local azimuthal angle of the umbral core magnetic fields for the 
same selection as in Figure]^ The azimuthal angles are defined as in Figure]^ 


where ANg corresponds the number of umbrae within a bin of inclination angles 
between 9q and 6i. The quantity A A/A gives the relative area of each inclination 
annulus given by the integral calculus for a unit sphere as 

AA/A = 27r(— cos 01 + cos0o)/47r = (— cos0i + cos0o)/2 (9) 

Finally, to normalize the curve by the area under the curve, we divide by p, 
which is the total sum of ANg/AA/A for all 9s. These distributions, given in 
Figure confirm the preference for umbral core magnetic fields to be radial; 
though, they also show the large spread of inclination angles about the radial 
direction. 

Figure]^ shows the histograms of the disambiguated magnetic azimuthal an¬ 
gles for the four categories. An azimuthal angle of 0° describes a vector directed 
towards solar west, as indicated in Figure]^ for positive and negative polarity 
spots with a diverging and converging magnetic field, respectively. While the 
later histograms exhibit a large spread, the umbral core magnetic fields of the 
trailing spots exhibit a slight preference to be located following Joy’s law. In the 
northern hemisphere, the trailing, positive-polarity spots have a mean azimuthal 
angle of —35.88°, consistent with a divergent field directed towards the South- 
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West. Likewise, the trailing, negative-polarity spots in southern hemisphere are 
consistent with a converging field directed away from the North-West direction. 
The leading polarities are not clearly located following Joy’s law. We speculate 
that these relationships are influenced by the asymmetric distribution of umbrae 
in the leading versus trailing polarities; that is, the leading polarity flux of a 
bipolar active region is typically more concentrated than the trailing polarity 
flmc (Fan, Fisher, and Deluca, 1993). 


6. Distributed Umbral Properties During The Solar Cycle 


The accurate characterization of the mean properties of sunspots requires a 
rigorous analysis of the distribution of these properties. Historically, this has 
not been a simple task due to a number of reasons: 1) the lack of instrumental 
stability and/or the variation of observing conditions, 2) selection bias intro¬ 
duced by manual sunspot identification and/or poor sampling, 3) low number 


statistics at solar minimum, and 4) the evolution of solar features. Bogdan 


et al. (1988) overcame the low number statistics by averaging over all solar 


cycles between 1917 and 1982 measured at Mount Wilson, and identified a log¬ 
normal distribution of umbral areas, whose average shape does not change over 
the solar cycle. Data from the Royal Freenwich Observatory support this log¬ 


normal distribution (Baumann and Solanki, 2005). Despite not being a full disk 


instrument, the Hinode/SP data benefits from its stable space-based platform, 
and our automated analysis eliminates user bias in selecting the umbrae. As 
long as our umbral size distribution compares well with historical measurements, 
we argue our distribution of core magnetic field strengths represents the most 
accurate determination to date. Previous distributions of umbral core magnetic 


field strengths have been presented by 

Livingston et al. 

(|2006 

), Rezaei, Beck, 

and Schmidt ( 

2012|), and |Livingston, Penn, and Svalgaard|( 

2012). 

We characterize the distribution of umbral property ‘x’ by its probability 


density function (PDF), (/(x), which for a continuous population is defined as 
the derivative of the cumulative distribution function (CDF), D{x) . Subject to 
the normalization requirement, the CDF states that 


r 

lim D{x) = lim / (j){x)dx = 1. 

X—^’OO X—>-00 / ^ 

J —oo 


( 10 ) 


Understanding the CDF as the probability that ‘x’ is less than some real-value 
within a population of size N, we can define the discretized CDF, D'{x), as 


lim D'{x) = lim 


1 AfV 

x^ao N Ax 

bins<ai 


Ax = 1, 


such that the discretized PDF takes the form: 

1 AN 
N Ax 


(/'(x) = 


( 11 ) 


( 12 ) 
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Figure 9. Discrete probability density functions (PDF) of umbral sizes (left panels) and core 
magnetic field strengths (right panels). The top plots correspond to the time period between 
November 2006 and December 2010, while the time period of the bottom plots is January 
2011 to November 2012. Continuous PDF fits are provided as indicated. See Table ^ for fit 
parameters for these time periods and for all the data. 


AN is the number of samples within a bin of width Ax centered on x. We 
determine 4>'{x) by counting the number of samples within a determined bin 
size and dividing that number by the bin size and the total number of samples 
in the population. This is similar to the umbral area spectrum definition given 
by Bogdan et al. (1988), except that we have introduced the normalization 
requirement. To compare our umbral area PDF with Bogdan’s et al., we divided 
their distribution functions by their scaling factor, /, and the number of samples 
in their umbra population. 

The discretized PDFs for umbral areas (given in units of microhemispheres) 
and core magnetic field strengths (in cgs units) in our Hinode/SP sample are pre¬ 
sented in Figure]^ Error bars representing the standard error {ax = y/AN/N/AX)^ 
are shown. The data are further split into two different time periods, November 
2006 to December 2010 and January 2011 to November 2012, which gives some 
indication of the distributions as a function of solar cycle amplitude. Using the 
regression analysis tool MPFIT (Markwardt, 2009), functional forms for the 
continuous PDFs are fit to the observed discrete PDFs and Icr parameter errors 
are derived. The goodness of fit is described by the normal reduced chi-squared 
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Table 1. Parameters for Umbral Magnetic Field Strength and Size Distributsion. Parameterl 
errors are shown between brackets in the second, third, and fourth columns. | 

Core Magnetic Field 

Log-Normal 

Time Period {Nspots) 
2006 - 2011 (179) 

2011 - 2013 (417) 

2006 - 2013 (596) 

<A> 

2062.64 (28.39) 
2040.97 (13.42) 
2065.67 (11.49) 

o-A 

1.0211 (0.0033) 
1.0161 (0.0012) 
1.0163 (0.0010) 

(l/Af)(AAr/Aa)^,,,, 
0.0012 (0.0001) 
0.0014 (0.0001) 
0.0014 (0.0001) 

0.78 

1.78 

2.21 

S ke w-Nor mal 





Time Period (Nspots) 

Location (^) 

Scale (cj) 

Skew (a) 

X^ 

2006 - 2011 (179) 

1830.97 (77.49) 

426.27 (64.1387) 

1.6229 (1.0056) 

0.82 

2011 - 2013 (417) 

1722.50 (18.24) 

501.47 (26.0172) 

5.0831 (0.9923) 

0.71 

2006 - 2013 (596) 

1730.31 (15.53) 

518.59 (22.1018) 

5.1874 (0.8081) 

0.87 

Umbral Areas 

Log-Normal 





Time Period {Nspots) 

<A> 

O-A 

(l/N){AN/Aa)roax 

X^ 

2006 - 2011 (179) 

0.0014 (0.0056) 

262.37 (716.11) 

17.813 (54.250) 

0.49 

2011 - 2013 (417) 

0.0011 (0.0032) 

170.74 (319.25) 

39.404 (97.371) 

1.45 

2006 - 2013 (596) 

0.0012 (0.0027) 

214.20 (320.80) 

27.300 (50.002) 

1.33 

Power Law 





Time Period (Nspots) 

Constant (C) 

Exponent (k) 

- 

X^ 

2006 - 2011 (179) 

0.4126 (0.0396) 

-1.5582 (0.0522) 

- 

0.58 

2011 - 2013 (417) 

0.4587 (0.0283) 

-1.6920 (0.0386) 

- 

1.49 

2006 - 2013 (596) 

0.4442 (0.0229) 

-1.6320 (0.0299) 

- 

1.46 


parameter, ■ To compare with Bogdan et al. (1988), a log-normal function is 


ascribed to the umbral area distribution, using a similar form, i.e. 


In 


1 AiV 
N Ax 


(In a — ln(a))^ 
2 In ax 


In(l^) 
N Ax J, 


(13) 


Additionally, we fit a power law {(j)'{x) = Cx^) to the observed umbral area 
PDFs. All function parameters and errors are given in Table [T| We find that 
the distribution of umbral sizes compares well with Bogdan et al. (1988) for 


both time periods, in addition to the full population. The log-normal description 
works well, with slightly better fits than using power laws for all time periods. 
We do not expect the distribution to match perfectly with that of|Bogdan et al. 


since Hinode probes a different portion of the solar continuum and our umbra- 
penumbral intensity definition is likely different. Yet, the log-normal shape of 
our distribution confirms that our selection is an adequate representation of the 
distribution of umbral parameters on the Sun. 

The PDF of the core magnetic field strength should be the transformation 
of the size PDF subject to the relationship between size and magnetic field 
discussed in Section]^ Previous investigations of the umbral area-magnetic field 
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relationship have found a linear relationship (Rezaei, Beck, and Schmidt, 20121; 


though, our high-resolution data displays a relationship that deviates somewhat 
from linearity. If the relationship is linear, one would expect a log-normal distri¬ 
bution in the core magnetic field strengths. Since our relationship is not linear, 
we seek an additional functional form to describe the magnetic field PDF. We 


introduce the skew-normal function described by Azzalini (1985): 


f(x,a) = 2(j){x)<^{ax), 


(14) 


where and $(a:) are the standard normal (Gaussian) probability density 
and cumulative distribution functions defined as 


Hx) = 


1 






$(c 


r 

= / 

j —( 


(j)(t)dt = - 


1 -b erf 


- C) 

V2 


UJ 


(15) 

(16) 


Through this definition, the normal parameters describing the Gaussian density 
function have different meanings. w, a are known as the location, scale, and 
shape (or skew) parameters, respectively. When a = 0, f{x,a = 0) is equivalent 
to a standard normal Gaussian distribution, where the ^ and ui give, respectively, 
the mean and the standard deviation. Although log-normal fits still fit well to 
the data (y^ < 2.21), skew-normal fits seem to better describe (y^ < 0.87) the 
core magnetic held PDF for all time periods (see Table . Furthermore, the 
skew parameter gives us an idea of a non-Gaussian behavior. The skewness of 
the distribution for the later years (a = 5.08) is higher than that of the earlier 
years (a = 1.62) in our sample; perhaps due in part to the combined effect of the 
nonlinear relationship between magnetic held strength and size and the greater 
raw number of larger spots in the years 2011 and 2012. Each ht supports a non¬ 
zero level of skewness in the PDF, as opposed to the results of |Rezaei, Beck, and| 
Schmidt (2012) and Livingston, Penn, and Svalgaard (2012). 


Skewness changes in the umbral core magnetic held PDF may imply variations 
in the mean magnetic held strength over time. While our time series is short, 
the skewness parameters indicate a higher skew when the level of activity is 
greater (i.e. when the number of sunspots is larger). This is consistent with 
the presumed steady size distribution of umbrae and the nonlinear relationship 
between size and magnetic held strength. When the cycle amplitude is larger, an 
identically shaped size distribution produces a greater number of larger spots, 
although the relative number of these spots does not change signihcantly. The 
increased number of larger spots reduces the level of sampling noise in the upper 
tail of the magnetic held strength distribution such that the non-linear relation¬ 
ship between umbral size and core magnetic held strength better inhuences the 
determined magnetic held strength distribution. The associated change in the 
arithmetic mean can be calculated by the expectation value for a skew-normal 
distribution given by 


E[x] = ^ + i 


vT" 


(17) 
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Figure 10. (top) Umbral core magnetic field strength of all 7530 identified umbrae (data 
points) as a function of time, (middle) The LOS component of the magnetic field divided 
by fi (i.e. the vertical (or radial) field approximation) for all spots, (bottom) The difference 
between LOS estimates for the total core magnetic field strength and the true core magnetic 
field strength. Blue squares denote the mean field strength for half-year bins, with error bars 
indicating the Icr standard error of the mean. 


The arithmetic mean for the three skew-normal functional fits to the magnetic 
field PDF give 2120.53, 2115.09, and 2136.60 G, respectively. Thus, while the 
magnetic field PDF shape changes, the mean magnetic field does not significantly 
vary between the earlier and the later years of our sample. 

Finally, we investigate the possibility of long-term secular trends in the mean 
umbral magnetic field strength, as reported by Penn and Livingston (2011). 
We are also interested in the utility of LOS magnetic flux measurements to 
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investigate trends in the mean magnetic field strength over the solar cycle, when 
one applies the vertical (or radial) field assumption. Figure 10 shows the core 
magnetic field strength and the core LOS approximation for the magnetic field 
strength, as a function of time, for all the identified umbrae. The difference 
of the two measurements is also plotted. A linear ht to the mean magnetic 
field strength binned in half-year bins as a function of the year after 2006 {i.e., 
{/{year — 2006) = m x {year — 2006) -I- b) results in a negligible gradient. The 
fitted gradient returns m = (1.023 ±2.13) G yr“^, with b = (1929.45 ± 11.39)G. 
Using the LOS values, the fitted line is described by m = (5.76 ± 3.33) G yr“^, 
b = (1827.84 ± 17.54) G. Gompare th is with the average gradient shown by 
Livingston, Penn, and Svalgaard (2012) of (—46 ± 6) G yr“^ between 1998 and 


2010. None of our measurements give a long- term decreasing trend. We do find 
short-term trends in the data implying possible solar cycle variations, which we 
claim can occur alongside a steady size distribution. For example,between July 
2009 and November 2012, the mean core magnetic field strength shows a gradient 
of (—24.99±5.39( G yr“^ while the mean ISlosI/m values display a gradient of 
(—39.55± 8.66) G yr“^. The LOS measurements clearly overestimate the trend, 
which is consistent with our discussion on the deviation of the magnetic field 
vector in dark umbral cores from the strictly radial approximation. We do not 
think that the trends shown between July 2009 and November 2012 mean or 
confirm a long-term decrease in the mean umbral magnetic held strength, i.e., 
on a time scale longer than a solar cycle. Rather, it suggests the mean magnetic 
held strength may vary over the solar cycle. We hnd no evidence for a long-term 
secular decrease. 


7. Discussion 

Our collection of umbrae observed with by the Hinode/SP ohers the highest 
resolution account of the thermal and magnetic properties of dark umbral cores 
to date. As studies continue to use the collective properties of umbrae to study 
the variations of the mean sunspot structure and formation over the solar cycle, 
it is critical to address the complexities involved in relating the various properties 
of sunspots to the variations of these properties over the solar cycle. Furthermore, 
understanding the limitations of the various data sources used is very important. 

First, and foremost, we claim that our selection of umbrae adequately char¬ 
acterizes the full spectrum of umbrae on the Sun, despite the fact that the 
Hinode/SP is not a full-disk instrument. We claim this because the determined 
PDF of umbral areas is consistent with previous studies, both for the whole 
portions of the solar cycle. Although |Nagovitsyn, Pevtsov, | 
) suggest that the distribution of the total sunspot areas 
(including penumbrae) varies with time, we find no support for changes in the 
umbral area PDF with time. The classification and identification of penumbra in 
a synoptic data set is quite complicated, as penumbrae associated with different 
umbrae are often merged and/or too close to each other. The identification of 
individual umbrae is more straight-forward and reliable. 

Both the normalized intensity and the magnetic held strength measured in 
the darkest region in each umbra show a clear correlation with the total size of 


cycle and for different 


and Livingston (2012 
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the umbra. The relationships are not linear. Consequently, the distribution of 
these parameters are not expected to follow the same distribution as the umbral 
areas. Furthermore, while the umbral area PDF remains constant, its amplitude 
varies as a function of solar cycle phase and as a function of the total solar cycle 
amplitude. This leads to a lower or greater number of large spots from cycle to 
cycle (with no change in the relative number of these spots). This, we expect, 
explains the variations in the maximum magnetic field strength of the umbrae 
over the solar cycle reported by Pevtsov et al. (2011). The size of the largest 
sunspots increases and, then, decreases as a function of the solar cycle as shown 
by Watson, Fletcher, and Marshall (2011), which in turn induces a solar cycle 
variation in the average umbral field strength inferred from the largest sunspots 
only. 

Due to the nonlinear relationship between magnetic field strength and umbral 
area, the magnetic field PDF may show small variations as a function of of the 
solar cycle evolution simply due to the change in the cycle amplitude. This fact, 
coupled with the fact that the distribution of umbral areas on the Sun at a 
particular time may deviate significantly from the average umbral area PDF, 
can lead to variations in the mean umbral magnetic field strength as a function 
of time. Such changes would not necessarily imply a significant deviation the 


solar dynamo activity, 

as suggested by 

Livingston, Penn, and Svalgaard 

(2012). 

Still, in opposition to 

Livingston, Penn, and Svalgaard ( 

2012 

), we do not find 


evidence for any long-term decrease in the mean magnetic field strength for the 
analyzed period. 

Furthermore, our determined magnetic field PDF does not support the normal 
distribution reported by Livingston, Penn, and Svalgaard (2012). We note that 
our work has a considerably larger sample of umbrae. According to |Livingston| 
et al. (2006) and Livingston, Penn, and Svalgaard (2012), the McMath-Pierce 


Telescope 1565 nm umbral magnetic field data set includes 900 sunspots, mea¬ 
sured across the entire solar disk between 1998 and 2005, and a total of 300 
observing days between 1998 and 2011, amounting to an average sampling rate 
of 6%. Meanwhile, our Hinode/SP data set includes a total of 7530 umbrae 
between 2006 and 2012 during a time of low solar cycle ampltiude; 596 of 
these 7530 umbrae are unique, as they are measured close to solar meridian. 


Livingston, Penn, and Svalgaard (2012) do not record size information on the 


measured sunspots that can be compared with the previous observations of the 
umbral size spectrum. Their observations of a observations of a core magnetic 
field distribution, whose shape remains unchanged while it mean decreases in 
time is inconsistent with our Hinode/SP measurements. While we assert that 
the magnetic field PDF can vary over the solar cycle, this variation is restricted 
primarily to the upper tail of the distribution. Due to the relative sparseness 
of the measured spots between 1998 and 2002 in the McMath-Pierce Telescope 
data set, the reported yearly change in the arithmetic mean is more suggestive 
of a sampling bias. 

To conclude, the collective study of simple measures of umbral properties 
provides constraints on the processes of continuum sunspot formation, and the 
connection of these measures with the bulk activity of the solar dynamo. Un¬ 
derstanding the variation of these mean properties over time is complicated 


SQLA: schad_revised2.tex; 22 May 2015; 0:16; p. 22 





































On the Collective Vector Magnetic Field Structure of Dark Umbral Cores 


because of the inherently nonlinear growth and evolution of individual sunspots. 
Furthermore, we must be careful to not assume that the properties of large 
umbrae to hold for smaller dark features, as in the example of the radial field 
approximation for dark umbral cores. Yet, viewing sunspot umbral cores as a 
continuum of sunspot formation seems to be valid as the relationship between 
magnetic field and temperature provide tight coherent relationships, albeit non¬ 
linear and reflective of involved radiative MHD processes, that might be helpful 
in comparison with advanced modeling efforts such as Rempel et ai (2009). 
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